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A  PRELIMINARY  WEATHER-CROP  YIELD  MODEL 
FOR  SOVIET  GRAIN 


Valentine  Zabijaka* 


Introduction 


Much  research  has  been  directed  toward  the  problem  of  relating  crop  yields 
to  weather.  The  value  of  a  model  that  would  show  the  effect  of  weather 
factors,  such  as  temperature  and  precipitation,  on  the  yield  of  a  specific 
crop  can  be  easily  justified.  Perhaps  the  most  important  use  of  a  predic¬ 
tive  weather-crop  yield  model  would  be  in  evaluating  the  potential  produc¬ 
tion  of  a  given  crop. 

The  analysis  that  follows  is  a  brief  review  of  preliminary  research  aimed 
at  estimating  the  effects  of  weather  (temperature,  precipitation,  and  soil 
moisture),  fertilizer  use,  and,  in  some  instances,  time  on  wheat  yields  in 
the  USSR.  Regression  analysis  was  used  to  determine  weather  variables  that 
significantly  influence  crop  yields.  The  resulting  regression  equations 
are  used  to  "predict"  wheat  yields  in  the  USSR  for  1972  and  1973.  Basic 
research  was  completed  in  August  1973.  Minor  revisions  include  incor¬ 
poration  of  data  on  plan  fulfillment,  published  in  January  1974. 

Procedures 


Data :  Weather  regions  of  the  USSR  (see  map)  were  selected  as  basic  geo¬ 
graphical  areas  for  computational  purposes.  The  27  weather  regions  do 
not  cover  the  USSR  but  do  include  about  90  percent  of  total  wheat  pro¬ 
duction.  Final  estimates  of  total  USSR  wheat  production  were  adjusted 
to  include  the  residual  production,  that  is,  the  10  percent  in  areas 
outside  the  weather  regions. 

Data  on  average  precipitation  and  temperature  were  available  for  each 
weather  region  and  could  be  used  as  given.  0ata  on  production  and  yields 
of  winter  and  spring  wheat  are  derived  from  official  Soviet  statistics 
for  crop  regions  (which  in  several  instances  correspond  to  weather  regions) 
and  oblasts.  Computations  were  based  on  the  dominant  type  of  wheat  grown 
in  the  given  weather  region.  In  a  few  regions,  such  as  those  in  the  Volga 
basin,  computations  were  performed  for  both  types.  Data  for  the  time 
period  1960-71  were  used  to  estimate  weather-yie Ld  relationships. 


*Based  on  a  paper  presented  by  Valentine  Zabijaka,  formerly  Economic 
Research  Service,  U.S.  Department  of  Agriculture,  at  a  Seminar  sponsored 
by  the  Foreign  Demand  and  Competition  Division,  Economic  Research  Service, 
September  14,  1973.  Mr.  Zabijaka  is  now  with  tha  U.S.  Department  of  Commerce. 
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Varlables:  The  following  independent  variables  were  used  in  constructing 
the  regression  equations: 

(1)  Precipitation  -  Total  precipitation  during  a  month 
in  millimeters. 

(2)  Soil  moisture  -  Estimated  by  a  simplified  Thornthwaite 
formula  for  the  last  day  of  each  month. 

(3)  Temperature  -  The  simple  monthly  average  of  daily  mean 
temperatures . 

(A)  Fertilizer  -  Statistical  data  is  available  on  total 
chemical  fertilizer  used  in  a  given  republic;  no  data 
is  available  on  fertilizer  application  to  individual 
crops, or  by  region. 

(5)  Time  -  In  a  few  regions  time,  in  years  was  used  as  a 
variable  in  place  of  fertilizer.  Fertilizer  use  and 
time  are  closely  correlated. 

Regression  techniques:  Step-wise  multiple  linear  regression  was  used,  with 
computer  selection  of  variables.  The  independent  variable  first  entered  in 
the  computation  was  selected  on  the  basis  of  highest  correlation  with  the 
dependent  variable  (yield) .  Subsequent  independent  variables  were  entered 
according  to  the  extent  to  which  they  explain  the  residual  fluctuations  in 
the  dependent  variable.  Only  those  variables  with  significant  t-values  are 
included  in  the  final  equation.  A  weather-yield  equation  was  directly  es¬ 
timated  for  each  weather  region. 

Results 


Wheat  yield  regression  equations  for  the  27  USSR  weather  regions  are  shown 
in  Table  1.  "Predictions"  of  Soviet  wheat  yields  in  each  region  for  1972 
and  1973  are  shown  in  Table  2.  Explanatory  notes  follow  each  table. 

Fit  of  the  equations  (Table  1)  appears  to  be  good.  Coefficients  of  deter¬ 
mination  (R^)  are  high;  21  are  .85  and  above.  Relatively  poor  fits  were 
obtained  for  two  regions  (South  Urals  and  West  Kazakhstan).  By  omitting 
1965,  a  poor  year,  much  better  fits  are  obtained.  Both  alternatives  are 
shown.  It  is  recognized  that  the  good  fit:  of  the  regression  equations 
may  be,  in  part,  a  result  of  the  limited  number  of  years  of  data  in  re¬ 
lation  to  the  number  of  independent  variables  used,  i.e.  a  small  number 
of  degrees  of  freedom. 

Fertilizer  appears  to  be  an  important  explanatory  factor  in  nearly  all  of 
the  winter  vheat  regions  but  in  only  one  or  two  of  the  spring  wheat  regions. 


. 
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A  time  trend  term  was  used  in  one  spring  wheat  region  (West  Siberia) .  Time 
trends  could  have  been  used  in  other  regions  instead  of  fertilizer  as  the 
two  factors  are  closely  correlated.  Fall  soil  moisture  (October  or  November) 
appears  to  be  the  most  important  weather-type  variable  in  the  winter  wheat 
equations,  while  June  conditions  (soil  moisture,  precipitation,  and  tem¬ 
perature)  explained  more  of  the  yield  variation  in  the  spring  wheat  regions. 

Using  the  regression  technique  indicated,  an  overall  USSR  winter  wheat  yield 
27.3  quintals  per  hectare*  (40 . 6  bushels  per  acre)  was  predicted  for  1973. 
This  can  be  compared  to  overall  USSR  winter  wheat  yields  (derived  from 
Soviet  statistics)  ranging  from  12.9  to  23.1  quintals  per  hectare  (19.2 
to  34.3  bushels  per  acre)  during  1960-71.  Using  regression  techniques, 
an  overall  yield  of  11.2  quintals  per  hectare  (16.7  bushels  per  acre) 
was  predicted  for  spring  wheat  in  1973.**  The  predicted  yield  of  spring 
wheat  is  within  the  range  of  overall  USSR  spring  wheat  yields,  which  in 
the  1960-71  period  ranged  from  5.5  to  12.2  quin  ;als  per  hectare  (8.2  to 
18.1  bushels  per  acre). 


Total  USSR  wheat  production  for  1973,  based  on  :he  previously  described 
regression  technique  and  incorporating  preliminary  and  revised  area  es¬ 
timates  is  shown  below. 


Preliminary 

Revised 

Area 

Prod . 

Area 

Yield 

Prod . 

Mil.  ha. 

Mil.  tons 

Mi  L .  ha . 

Quintals/ha. 

Mil.  tons 

Winter  wheat 

16.0 

43.7 

L8.3 

26.2 

48.0 

Spring  wheat 

46.5 

52.0 

44.7 

13.8 

61.7 

17  Official 

62.5 

Soviet  figur 

95.7 

e . 

63.0 

17.4 

109.7  1/ 

Research  is  continuing  and  efforts  will  b<  made  to  refine  the  regression 
equations  through  input  of  new  data  and  u;  e  of  other  regression  and 
statistical  techniques. 


*Conversion  equivalents:  One  metric  quintal  (centner)  =  220.46  pounds 

One  hectare  =  2.471  acres 

To  convert  metric  quintals  (centners)  of  wheat  per  hectare  to  bushels  per 
acre,  multiply  by  1.487. 

**A  yield  of  9.7  quintals  per  hectare  would  result  if  alternative 
equations  omitting  1965  data  were  used  for  the  South  Urals  and  West 
Kazakhstan  regions. 
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56068  1?  60 


Weather  Regions 


Number  Name 


1 

Baltic 

2 

Belorussia 

3 

West  Ukraine 

4 

North  Central  Ukraine 

5 

Northeast  Ukraine 

6 

East  Ukraine 

7 

South  Ukraine 

8 

Moldavia 

9 

Krasnodar 

10 

Northeast  Caucasus 

11 

West  Blacksoil 

12 

East  Blacksoil 

13 

Central 

14 

Volga-Vyatsk 

15 

Upper  Volga 

16 

Middle  Volga 

17 

Lower  Volga 

18 

Northwest  Urals 

19 

South  Urals 

20 

Northeast  Urals 

21 

West  Kazakhstan 

22 

Kustanai 

23 

Tselinograd 

24 

North  Kazakhstan 

25 

Pavlodar 

26 

West  Siberia 

27 

Altai  Krai 

WHEAT  YIELD  REGRESSION  RUNS  FOR  THE  USSR  WEATHER  REGIONS 
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table  l  WHEAT  YIELD  REGRESSION  RUNT  OR  THE  USSR  WEATHER  REGIONS 


Assumed,  yield.  No  yield  data  available  by  type  of  wheat 

Regression  equation  without  1965  data. 

lata  for  selected  years:  i960, 1961, 1965-68,  and  1971* 


WHEAT  YIELD  REGRESSION  hoi-TS  FOR  THE  USSR  WEATHER  REGIONS 
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Table  1:  Explanatory  Notes 


Column 

Remarks 

1  USSR  Weather  Regions  :  As  shown  on  map.  Regions  delineate  areas 

with  similar  soil,  crop,  and  weather  characteristics.  The  percentage 
figure  indicates  the  relative  importance  of  the  region  in  terms  of  its 
wheat  production  (1971) . 

2.  "Wheat  type'':  Winter  wheat  (WW)  and  spring  wheat  (SW).  The  dominant 
type  of  wheat  grown  in  each  region  is  noted.  (Both  types  are  shown  for 
some  regions.) 

3.  "Regression  equation":  The  yield  regression  equations  are  expressed 
in  quintals  per  hectare.  Independent  variables  are  shown  in  the  order 
of  their  entry  into  the  regression  equation.  The  following  independent 
variables  were  used: 

a.  Temperature  ("T",  followed  by  the  abbreviation  of  the  month^-- 
The  simple  monthly  average  of  daily  mean  temperatures  (Celsiils, 
to  one  decimal  place) . 

b.  Precipitation  ("P",  followed  by  the  abbreviation  of  the  month) — 

The  total  precipitation  during  a  month  in  millimeters. 

c.  Soil  moisture  ("S",  followed  by  the  abbreviation  of  the  month) — 
Estimated  by  a  simplified  Thornthwaite  formula  for  the  last  day 
of  the  month. 

d.  Fertilizer  ("FT") — Based  on  the  estimated  amount  of  fertilizer  used 
per  hectare  of  cultivated  land  (kilograms  per  hectare). 

e.  Time  ("TIME") — Annual  trend  used  in  place  of  fertilizer  in  a  few 
regions.  Fertilizer  use  and  time  are  closely  correlated. 

4.  "r2":  Coefficient  of  determination.  Expresses  the  percent  of  total 
fluctuations  in  wheat  yield  explained  by  the  regression  equation. 

5.  "S.E."  [S.E.E.] :  Standard  error  of  estimate.  Measures  the  accuracy 
of  estimates  obtained  through  the  use  of  the  regression  equation  (68% 
of  all  cases  lie  within  ^  one  standard  error  of  estimate;  95%  within 
-2  S.E.E.;  99.7%  within  ±3  S.E.E.).  The  S.E.E.  refers  to  estimates 
within  the  time  period  of  the  regression;  the  expected  error  of  fore¬ 
cast  yields  would  be  greater. 

6.  "Mean":  Observed  average  wheat  yield  in  the  given  weather  region 
during  the  study  period  (usually  1960-71). 


. 
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7.  "S.D.":  Standard  deviation,  a  measure  of  variation  around 

the  observed  mean,  similar  to  the  standard  error  of  estimate 
(68%  of  all  cases  lie  within  a  range  of  -  one  standard  devi¬ 
ation,  and  so  on). 
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Table  2 :  Explanatory  Notes 


Column 

Remarks 

1  "Weather  region":  As  shown  on  map.  Regie  ns  delineate  areas 
with  similar  soil,  crop,  and  weather  characteristics. 

2  "Type  of  wheat":  Winter  wheat  (WW)  and  spring  wheat  (SW) . 

The  dominant  type  of  wheat  grown  in  each  region  is  noted. 

(Both  types  are  shown  for  some  regions.) 

3,4  "Projection  1972,  1973":  Yield  projections  for  1972  and 
1973  in  quintals  per  hectare. 

5  "r2":  Coefficient  of  determination*  Expresses  the  percent  of 
total  fluctuations  explained  by  the  regression  equation. 

6  "Standard  error":  Standard  error  of  estimite  (S.E.E.).  Mea¬ 
sures  the  accuracy  of  estimates  obtained  through  the  use  of 
the  regression  equation  (68%  of  all  cases  lie  within  -  one 
standard  error  of  estimate;  95%  within  i 2  S.E.E.;97.7%  within 
±3  S.E.E.).  The  S.E.E.  refers  to  estimates  within  the  time 
period  of  the  regression;  the  expected  error  of  forecast  yields 
would  be  greater. 

7,8  "1960-1971  range":  Wheat  yield  range  during  the  study  period 

(usually  1960-71)  in  quintals  per  hectare. 

9  "Average  yield":  Average  (mean)  wheat  yield  over  the  study 

period  (usually  1960-71)  in  quintals  per  hectare. 

10  "Standard  deviation":  A  measure  of  variation  around  the  average 
(mean)  wheat  yield,  similar  to  the  standard  error  of  estimate 
(68%  of  all  cases  lie  within  a  range  of  ^  one  standard  deviation, 
and  so  on) . 

11  "Comments":  Record  yield  projections  noted. 
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